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Abstract 

Solidification  modeling  has  been  utilized  to  predict  local  thermal 
conditions  during  solidification  for  a  broad  range  of  geometrical 
configurations,  alloy  compositions  and  heat-extraction  conditions. 
Three  different  mold  configurations  were  evaluated  for  three  alloy 
compositions  using  both  conventional  and  high-gradient 
directional  solidification  processes.  The  high-gradient  directional 
solidification  process  investigated  was  the  Liquid  Metal  Cooling 
(LMC)  process  that  utilizes  a  liquid-metal  coolant  in  the  cold  zone 
of  the  directional- solidification  furnace.  Predictions  of 
solidification  conditions  have  been  compared  to  classical  defect 
formation  criteria  to  determine  solidification-model  applicability. 
The  classical  defect  maps  and  dendrite-growth  models  have  been 
extended  to  accommodate  curvature  of  the  solid-liquid  interface, 
which  has  not  been  well  understood  to  date.  Experiments  have 
been  conducted  to  validate  model  predictions  and  improve  the 
understanding  of  the  role  of  solid-liquid  interface  curvature  on 
dendrite-growth  morphology.  In  particular,  lateral  growth  - 
secondary  dendrite  arms  overgrowing  primary  arms  well -aligned 
with  the  withdrawal  direction  -  has  been  identified  as  a  precursor 
to  misoriented  grain  formation.  A  means  to  predict  preferred 
process  conditions  for  directional  solidification  has  been 
developed  using  the  axial  thermal  gradient  at  the  surface  of  the 
casting  during  solidification.  The  utilization  of  the  surface 
maximum  axial  thermal  gradient  as  a  means  to  identify  preferred 
processing  conditions  is  applicable  to  a  range  of  solidification 
conditions,  and  accounts  for  changes  in  casting  geometry,  alloy 
composition  or  degree  of  heat-extraction.  This  optimization 
technique  has  been  demonstrated  for  an  atypical  casting 
configuration  and  applied  to  a  complex  geometry,  in  which  time- 
dependent  process  conditions  were  required  to  maintain  desired 
single-crystal  growth. 


Introduction 

Single-crystal  nickel-base  superalloys  have  been  utilized  for  high- 
pressure  turbine  blades  for  decades  due  to  their  exceptional  high- 
temperature  performance.  These  materials  are  currently  produced 
by  directional  solidification  via  the  Bridgman  process  (Figure  la). 
During  this  process,  an  investment  casting  is  withdrawn  from  the 
hot  section  to  the  cold  section  of  a  vacuum  furnace  [1-2].  The 
process  relies  on  radiation  heating  in  the  hot  zone  and  radiation 
cooling  in  the  cold  zone  to  maintain  a  relatively  flat  solid-liquid 
interface  that  proceeds  along  the  casting  as  it  is  withdrawn  in  the 
furnace.  For  large  or  complex  geometries  it  is  difficult  to 
maintain  a  planar-growth  front  and  a  small  processing  window 
exists  due  to  breakdown  in  single-crystal  growth  from  inefficient 
radiation  cooling. 

The  Liquid  Metal  Cooling  (LMC)  directional -solidification 
process  has  been  developed  and  addresses  the  need  for  enhanced 
heat  extraction  (Figure  lb).  The  process  utilizes  a  liquid-metal 


coolant  in  the  cold- section  of  the  furnace  to  efficiently  extract  heat 
from  the  casting  [3,4].  The  process  has  demonstrated  3-9x 
increase  in  cooling  rate,  corresponding  to  a  50  pet  reduction  in 
micro  structure  scale  with  reduced  defect  occurrence  [4-18].  The 
refined  micro  structure  has  resulted  in  a  3-7x  improvement  in 
fatigue  life  on  simple  geometries  [18].  However,  the  enhanced 
lateral  heat  extraction  from  the  surface  of  the  mold  can  promote 
solidification-front  curvature,  which  promotes  the  formation  of 
misoriented  grains  [16].  The  effect  of  solid-liquid  interface 
curvature  on  dendrite  morphology  and  on  the  propensity  to  form 
defects  is  not  well  understood. 


(b) 


Figure  1 :  Schematic  of  the  (a)  Bridgman  and  (b)  LMC  processes 

[8]. 


Previous  LMC  experiments  have  observed  the  formation  of  a 
dendrite  morphology  referred  to  as  lateral  growth  (Figure  2)  [16]. 
It  was  demonstrated  that  lateral  growth  occurred  when  the  local 
thermal  gradient  was  better  aligned  with  a  secondary  dendrite  arm 
than  the  primary  dendrite  arm,  causing  the  secondary  arm  to  grow 
faster  than  and  overgrow  primary  dendrite  arms  well-aligned  with 
the  withdrawal  direction  [16].  The  lateral  growth  regimes  were 
large,  in  some  cases  completely  growing  to  the  center  of  bars  [16]. 
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Figure  2:  Micrographs  of  lateral  growth  (a)  transverse  and  (b) 
longitudinal  to  the  growth  direction  [16]. 


Several  previous  experimental  investigations  have  been  conducted 
to  develop  relationships  between  macroscopic  thermal  conditions 
and  dendrite  morphology  [19-23].  The  dendrite  morphology  has 
been  associated  with  the  thermal  gradient  during  solidification  and 
the  solidification  rate  (Figure  3).  Regimes  of  axial  growth, 
freckle  formation,  nucleated  grains  and  complete  breakdown  have 
been  defined.  The  underlying  mechanisms  that  control 
microstructure  scale  have  also  been  defined  in  similar  analyses 
that  demonstrate  refinement  with  increased  thermal  gradient  and 
solidification  rate.  Due  to  limitations  of  the  experimental 
techniques,  an  assumption  has  been  made  that  the  solid-liquid 
interface  is  flat  [19-21].  Thus,  classical  defect  maps  do  not 
currently  account  for  curvature  of  the  solid-liquid  interface.  In 
addition,  the  dendrite-growth  models  are  incapable  of 
incorporating  curvature  into  predictions  of  micro  structure  scale. 
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Figure  3:  Adapted  map  of  classical  defect-formation  criteria  with 
superimposed  schematics  of  dendrite  morphology  for  directional 
solidification  [19-20]. 

Continuum  solidification  models  have  been  used  to  develop  trends 
in  directional- solidification  behavior  based  on  broad  changes  in 
processing  conditions  [12-17,24-25].  The  ability  for  process 
simulations  to  accurately  predict  changes  in  local  thermal 


conditions  during  solidification  due  to  subtle  variations  in  process 
conditions  has  only  recently  been  possible,  especially  for  relevant 
geometrical  configurations.  State-of-the-art  solidification  models 
are  capable  of  utilizing  time-,  temperature-  and  location- 
dependent  material  properties  and  boundary  conditions  that  enable 
accurate  predictions  of  solidification  behavior.  However, 
continued  advances  predicting  stress  and  incorporating  pressure 
predictions  into  boundary  conditions  are  still  desired. 

The  characterization  of  local  thermal  conditions  as  a  function  of 
withdrawal  rate  and  furnace  temperature  for  directional 
solidification  is  challenging  and  typically  requires  multiple 
casting  trials.  The  process  parameters  are  generally  dependent  on 
alloy  composition  and  casting  geometry,  complicating  the 
determination  of  preferred  processing  conditions.  The  intent  of 
this  research  is  to  utilize  fundamental  solidification  theory  and 
solidification  modeling  to  develop  a  geometry-independent  metric 
that  can  be  applied  to  specific  casting  geometries  and  alloys  to 
identify  preferred  processing  conditions  for  a  broad  range  of 
solidification  conditions  that  include  curvature  of  the  solid-liquid 
interface.  Detailed  casting  experiments  have  been  employed  to 
validate  solidification-modeling  predictions. 

Procedure 


Experimental  Methods 

Casting  trials  using  an  ALD  directional- solidification  furnace 
were  utilized  to  support  model  validation  and  verification. 
Additional  details  of  the  furnace  have  been  documented  elsewhere 
[16].  A  unique  aspect  of  this  furnace  is  that  it  is  capable  of 
directional  solidification  in  Bridgman  and  LMC  modes,  providing 
a  broad  range  of  heat  extraction  conditions.  A  key  aspect  of  the 
LMC  furnace  is  the  utilization  of  a  floating  ceramic  baffle  atop 
the  coolant,  in  this  case  tin,  in  order  to  shield  the  coolant  from 
radiation  heating  via  the  mold  heater.  The  floating  baffle  was  1.5 
cm  thick  and  comprised  of  mullite  beads  3.2  mm  in  diameter. 

The  molds  for  each  casting  configuration  were  supplied  by  PCC 
Airfoils.  A  seed  and  grain  selector  were  utilized  for  all  castings  to 
ensure  single-crystal  growth  of  the  proper  crystallographic 
orientation.  Type  C  thermocouples  were  utilized  to  profile  the  hot 
zone  of  the  furnace,  obtain  cooling  curves  at  the  mold-metal 
interface  and  validate  solidification  predictions.  A  hole  was 
drilled  through  the  mold  for  proper  positioning  of  the 
thermocouples  at  the  mold-metal  interface.  The  nominal  mold 
thickness  was  0.8  cm.  Three  commercial  single-crystal 
superalloys  were  utilized:  Rene  N4,  Rene  N5  and  CMSX-486. 
The  alloy  compositions  are  summarized  in  Table  1. 


Table  1:  Compositions  of  investigated  alloys  (in  wt.  pet.). 


Alloy 

Cr 

Co 

Mo 

W 

Ta 

Re 

A1 

Ti 

Hf 

C 

Other 

Ni 

Rene  N4 

9.8 

7.5 

1.5 

6.0 

4.8 

4.2 

3.5 

0.15 

0.05 

0.5Nb 

bal 

Rene  N5 

7.0 

7.5 

1.5 

5.0 

6.5 

3.0 

6.2 

0.15 

0.05 

0.01Y 

bal 

CMSX-486 

5.0 

9.3 

0.7 

8.6 

4.5 

3.0 

5.7 

0.7 

1.2 

0.07  0.015B 

bal 

Three  geometrical  configurations  were  investigated:  (a)  a  5 -bar 
cluster-mold  cast  of  Rene  N5  [17],  (b)  a  specially  designed  mold 
to  investigate  lateral  dendrite  growth  using  alloys  CMSX-486  and 
Rene  N4  [16]  and  (c)  a  tri-crystal  ring  segment  cast  using  alloy 
CMSX-486  (Figure  4)  [14-16]. 
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Figure  4:  Casting  geometries:  (a)  rotationally- symmetric  section 
of  a  5-bar  mold,  (b)  rotationally- symmetric  section  of  the  lateral- 
growth  mold,  (c)  tri-crystal  and  (d)  airfoil. 


The  5 -bar  cluster  consisted  of  5  bars  evenly  distributed  around  a 
support  post  with  the  center  of  the  bar  positioned  at  a  radius  of  5.1 
cm  from  the  center  of  the  cluster.  Each  bar  was  15.9  cm  long  and 
1.6  cm  in  diameter.  The  alloy  utilized  for  this  investigation  was 
Rene  N5.  The  withdrawal  rates  investigated  were  3.4  mm- min'1 
for  the  Bridgman  process  and  8.5,  12.7  and  21.2  mm-min"1  for  the 
LMC  process. 

The  lateral-growth  casting  utilized  a  5 -bar  cluster  previously 
described  with  a  modification  to  increase  the  relative  lateral  heat 
extraction.  Two  patches  of  woven  alumina  fabric,  herein  referred 
to  as  “wrap”,  were  cut  to  a  width  of  2.5  cm,  wetted  with  alumina 
cement  and  affixed  to  the  exterior  of  the  mold  encircling  the  bar 
and  centered  at  4.5  and  10.8  cm  above  the  bottom  of  the  bar. 

Each  bar  was  modified  with  a  different  wrap  thickness  of  0.5,  1.0 
and  1.4  cm,  including  one  bar  with  no  additional  wrap  and  one 
replicate  bar  with  a  wrap  thickness  of  1.0  cm.  The  additional 
wrap  was  intended  to  increase  the  lateral  heat  extraction  just 
above  the  wrap  relative  to  within  the  wrap  in  order  to  develop 
varying  degrees  of  solidification-front  inclination  at  the  surface  of 
the  bar  just  above  the  wrap.  Only  the  LMC  process  was  evaluated 
for  this  mold  configuration.  Withdrawal  rates  of  12.7,  16.9  and 
21.2  mm-min'1  were  investigated  for  alloy  CMSX-486,  while 
withdrawal  rates  of  8.5,  12.7  and  16.9  mm-min'1  were 
investigated  for  Rene  N4. 

The  tri-crystal  ring  segment  was  a  30-degree  segment  of  a  full 
ring  with  a  height  of  2.5  cm.  The  ring  segment  was  divided  into 
three  10-degree  segments  in  which  each  10-degree  segment  had  a 
seed  and  grain  selector  from  which  to  grow  a  different 
crystallographic  orientation.  The  primary  orientation  was  aligned 
with  the  axial  orientation  and  secondary  orientation  with  the  radial 
orientation  of  the  ring  segment.  This  configuration  resulted  in  10- 
degree  crystallographic  boundaries  between  the  three  segmented 
single-crystals.  Two  ring  thicknesses  (inner-  to  outer-diameter 
distance)  were  considered,  designated  as  thick  (5.1  cm)  and  thin 
(1.9  cm).  The  intent  of  this  investigation  was  to  assess  the 
feasibility  of  casting  multiply- seeded  ‘single-crystal’  rings  of 
different  thicknesses  using  conventional  and  high-gradient 
directional- solidification  processes.  Each  ring  thicknesses  was 
investigated  for  each  of  the  Bridgman  and  LMC  processes  with  a 
single  casting  trial  based  on  process  conditions  determined  using 
solidification  modeling  as  will  be  described  in  the  following 
section. 

Castings  were  sectioned  and  prepared  for  metallographic  analysis. 
The  Primary  Dendrite  Arm  Spacing  (PDAS)  and  Secondary 
Dendrite  Arm  Spacing  (SDAS)  were  measured  according  to  well- 
accepted  techniques  as  described  elsewhere  [16].  Relative 


differences  in  PDAS  and  SDAS  were  compared  with  relative 
changes  in  the  predicted  thermal  conditions  of  the  castings. 

Simulation  Methods 

A  finite-element  solidification  model,  ProCAST™  (a  trademark 
of  ESI  Group,  Inc.),  was  utilized  to  predict  spatial-dependent 
thermal  conditions  during  directional  solidification.  Details 
regarding  model  parameters  and  boundary  conditions  have  been 
summarized  in  Table  2  and  fully  described  elsewhere  [16].  The 
solidification  model  has  been  validated  from  four  different 
process  and  mold-configuration  conditions  encompassing  a  broad 
range  of  heat-extraction  conditions  by  comparison  of 
thermocouple  measurements  at  the  mold-metal  interface,  exterior 
mold  surface,  and  specific  locations  within  the  hot,  cold  and 
transition  zones  of  the  directional  solidification  furnace  [15-16]. 

Table  2:  Model  parameters  for  FE-based  solidification  modeling 


Thermophysical  Properties  Interface  Heat  Transfer  Coefficients 


Baffle  Therm.  Cond. 

1.0'  W/mK 

Sup  eralloy-C  eramic 

*750  W/m‘K 

Shell  Therm.  Cond. 

2.5  W/mK 

C  eramic-C  eramic 

500  W/mJK 

Shell  Thickness 

6  A  mm 

Tin-Ceramic 

4000  W/oTK 

Boundary7  Conditions 

Run  Parameters 

Shell  Emissivity 

0.4 

Maximum  Time  Step 

** 

Baffle  Emissivity 

0.2 

Spatial  Step 

0.5  mm 

Tin  Temperature 

250  °C 

Critical  Fraction  Solid 

0.4 

*  Initial  value  reported  exponential  decay  to  IOOW.Wk  during  solidification 
** Scaled  to  provide  equal  withdrawal  distance  (0.4mm)  between  steps 


A  key  boundary  condition  for  the  LMC  process  simulations  was 
location-dependent  heat-transfer  conditions  at  the  surface  of  the 
mold  that  varied  whether  the  particular  node  was  in  the  mold 
heater,  floating  baffle  or  liquid  coolant  regions  of  the  furnace 
(Table  3).  Initially,  the  entire  mold  was  in  the  mold-heater  region 
and  radiation  heating  occurred  from  the  heating  elements  through 
application  of  an  emissivity  on  the  surface  of  the  mold.  During 
withdrawal,  each  mold-surface  node  passed  through  the  floating 
baffle,  during  which  the  emissivity  was  turned  off  and  ceramic- 
ceramic  interface  heat  transfer  coefficient  applied.  The  ambient 
temperature  was  assigned  according  to  the  thermal  gradient 
occurring  through  the  thickness  of  the  floating  baffle.  Finally,  as 
each  node  passed  from  the  floating  baffle  into  the  coolant,  the 
interface  heat  transfer  coefficient  was  increased  to  a  value 
consistent  with  a  liquid  metal-ceramic  interface  and  the  ambient 
temperature  was  set  to  the  experimentally-observed  temperature 
of  the  tin  bath. 


Table  3:  Location-dependent  boundary  conditions  to  simulate 
interactions  of  floating  baffle  and  coolant  with  mold. 


Location 

Emissivity 

Interface  Heat 

Transfer  Coefficient 

(W-m2-^'1) 

Ambient 

Temperature  (°C) 

Above  Baffle 

0.4 

0 

— 

Within  Baffle 

0 

500 

1100-  325  °C 

Below  Baffle 

0 

4000 

250 

A  simulated  parametric  analysis  was  conducted  for  each  casting 
configuration  to  evaluate  local  thermal  conditions  as  a  function  of 
withdrawal  rate  conditions  (Figure  5  and  Table  4).  The  sensitivity 


3 

Approved  for  public  release;  distribution  unlimited 


of  thermal  conditions  to  bar  thickness,  mold-heater  temperature 
and  alloy  composition  was  also  evaluated  for  the  5 -bar  cluster 
configuration  (Figure  5).  Overall,  200  simulations  for  the  5 -bar 
cluster,  24  simulations  for  the  lateral-growth  cluster  and  20 
simulations  for  the  tri-crystal  configuration  were  analyzed. 
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Figure  5:  Simulated  two-parameter  analyses  for  5 -bar  cluster, 
withdrawal  rate  and  (a)  bar  thickness,  (B)  mold-heater 
temperature  and  (c)  alloy. 


Table  4:  Simulated  withdrawal  rates  for  the  lateral- growth  and  tri¬ 
crystal  ring-segment  casting  configuration. 


Mold  Configuration 

Alloy 

Process 

Withdrawal  Rates  (mm-min1) 

lateral  growth 

CMSX-486 

Rene  N4 

LMC 

LMC 

12.7,  16.9,21.2 

8.5,  12.7,  16.9 

thick  tri-  crystal 

CMSX-486 

Bridgman 

LMC 

0.8,  1.7,  2.5,  3.4,  4.2 

4.2,  5.9,  6.8,  7.6,  9.3 

thin  tri- crystal 

CMSX-486 

Bridgman 

LMC 

1.7,  2.5,  3.4,  4.2,  5.1 

5.1,  6.8,  8.5,  10.2,  11.9 

The  thermal  conditions  examined  include  (1)  the  thermal  gradient 
in  the  withdrawal  direction  (herein  referred  to  as  the  axial  thermal 
gradient),  (2)  the  solidification  rate  in  the  withdrawal  direction 
(herein  referred  to  as  the  axial  solidification  rate,  (3)  the  cooling 
rate  from  liquidus  to  solidus,  (4)  the  inclination  angle  of  the  solid- 
liquid-interface  (relative  to  horizontal),  (5)  the  solid-liquid- 
interface  position  relative  to  the  transition  region  in  the  furnace 
(the  top  of  the  fixed  baffle  in  Bridgman  and  top  of  the  floating 
baffle  in  the  LMC  process,  herein  referred  to  as  the  relative 
position,  Figure  6).  One  node  at  the  surface  and  center  of  each 
casting  configuration  was  evaluated:  at  the  mid-height  of  the  5 -bar 
cluster,  just  above  the  wrap  of  the  lateral-growth  casting  and  at  the 
mid-height  of  the  tri-crystal  ring  segment.  The  values  of  the 
thermal  conditions  were  analyzed  with  changes  in  process 
conditions. 


Figure  6:  Schematic  of  fraction- solid  contours  depicting  the  solid- 
liquid  interface  position  relative  to  the  transition  region  in  the 
furnace  for  the  (a)  Bridgman  and  (b)  LMC  processes. 

The  PDAS  (2 i)  and  SDAS  (22)  were  predicted  from  the  local 
thermal  conditions  according  to  well-accepted  dendrite-growth 
models,  Equations  1  and  2,  respectively  [20,21]. 


•N 

1 

CN 

II 

(i) 

^2  ~  ^2  (^)  1^3 

(2) 

In  these  equations,  Aj  and  A2  are  coefficients  dependent  on  alloy 
composition,  Gz  is  the  axial  thermal  gradient,  V  is  the 
solidification  rate  and  L  is  the  cooling  rate.  The  original 
development  of  the  dendrite  growth  model  for  primary  dendrite 
arm  spacing  assumed  that  the  thermal  gradient  was  normal  to  the 
growth  direction  and  utilized  the  thermal  gradient  and  alloy 
melting  range  to  approximate  the  length  of  the  parabolic  dendrite. 
Thus,  when  extending  the  model  to  curved  solid-liquid  interfaces, 
the  thermal  gradient  parallel  to  the  dendrite  growth  direction  -  the 
axial  thermal  gradient  -  should  only  be  considered.  Measured 
PDAS  and  SDAS  were  compared  to  predicted  PDAS  and  SDAS 
based  on  local  thermal  conditions. 

From  the  thermal-condition  data,  trends  of  thermal  conditions 
were  used  to  develop  a  means  to  quantify  the  quality  of 
directional  solidification.  Iteration  of  process  conditions  was  then 
employed  to  define  a  preferred  process  condition  for  each  casting 
configuration.  The  optimization  approached  maximized  heat 
extraction  while  also  limiting  the  probability  for  the  formation  of 
solidification  defects  for  a  given  casting  geometry  and  alloy.  The 
range  of  alloy  compositions  and  geometrical  configurations 
allowed  for  generalization  of  the  optimization  approach.  The 
thermal  conditions  were  also  compared  to  classical  solidification 
maps  and  defect  criteria  for  each  casting  configuration  to  assess 
the  validity  of  the  data  from  experience-based  processing  maps 
[19-20]. 

Finally,  the  optimization  protocol  was  applied  to  a  complex  airfoil 
geometry  (Figure  4).  Directional  solidification  of  a  4-airfoil 
cluster  of  CFM56-3  airfoils  was  simulated  for  the  LMC  process. 
A  range  of  constant  withdrawal  rates  were  evaluated:  6.4,  8.5, 
10.6,  12.7,  14.8,  16.9,  19.1  and  23.3  mm-min'1.  Thermal 
conditions  from  nodes  throughout  the  casting  were  analyzed  with 
respect  to  the  axial  position  of  the  node  along  the  blade  axis.  A 
technique  to  define  a  time-dependent  withdrawal  rate  from  a 
simulated  parametric  analysis  was  demonstrated. 
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Results 


Case  A:  five-bar  cluster 


Predictions  of  local  thermal  conditions  at  the  surface  of  a  bar  at  its 
mid-height  from  a  simulated  parametric  analysis  of  the  directional 
solidification  process  indicated  that  a  maximum  axial  thermal 
gradient  was  achieved  at  an  intermediate  withdrawal  rate  (Figure 
7).  The  axial  thermal  gradient  is  an  important  thermal  metric 
during  directional  solidification  because  the  primary  dendrite  arm 
spacing  is  controlled  by  the  thermal  gradient  in  the  growth  (axial) 
direction.  As  the  withdrawal  rate  increased,  the  solidification 
front  lowered  further  from  the  hot  zone  of  the  furnace.  At  a 
withdrawal  rate  of  11.9  mm-min"1,  the  solidification  front  was  0.2 
cm  above  the  floating  baffle,  in  the  transition  region  between  the 
hot  and  cold  zones  of  the  furnace,  and  a  maximum  axial  thermal 
gradient  of  82  °C-cm'1  was  achieved  at  the  surface  of  the  casting 
(Figure  7).  Any  further  increase  in  the  withdrawal  rate  caused 
significant  curvature  of  the  solid-liquid  interface,  thus  reducing 
the  axial  thermal  gradient  and  increasing  the  propensity  to  form 
defects  such  as  stray  grains  (Figure  7).  However,  the  axial 
thermal  gradient  in  the  center  of  the  casting  further  increased, 
causing  variations  in  the  solidification  conditions  from  the  surface 
to  the  center  of  the  casting.  The  cooling  rate  increased  with 
increasing  withdrawal  rate  in  all  cases,  despite  the  severely 
inclined  solid-liquid  interface  (Figure  7). 
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Figure  7 :  Predicted  sensitivity  of  thermal  conditions  to  withdrawal 
rate  at  the  surface  and  center  of  the  bar  at  its  mid-height  for  a  5- 
bar  cluster  composed  of  CMSX-486  processed  via  Bridgman  and 
LMC,  (a)  axial  thermal  gradient,  (b)  cooling  rate  and  (c) 
inclination  angle. 


Classical  defect-formation  criteria  have  been  extended  to  non- 
axial  heat  extraction  during  solidification  by  utilization  of  the 
axial  thermal  gradient  and  axial  solidification  rate.  Comparison 
between  experimental  observations  and  defect  criteria  validated 
the  use  of  the  axial  thermal  conditions,  as  will  be  discussed  in  the 
following  section.  The  predicted  axial  thermal  gradient  and  axial 
solidification  rate  of  each  of  the  simulated  castings  were 


superimposed  on  the  modified  classical  solidification  maps 
(Figure  8).  The  casting  conditions  that  achieved  the  maximum 
axial  thermal  gradient  were  within  the  directional  growth  zone  of 
the  solidification  maps.  An  increase  in  withdrawal  rate  from  this 
position  approached  the  thermal  conditions  consistent  with  stray 
grains,  in  which  a  change  in  dendrite  morphology  from  axial 
growth  to  lateral  growth  -  secondary  dendrite  arms  over-growing 
well-aligned  primary  arms  -  and  ultimately  stray  grain  formation 
was  observed.  Conversely,  a  decrease  in  the  withdrawal  rate 
approached  the  thermal  conditions  consistent  with  freckle 
formation. 
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Figure  8:  Defect  map  for  directional  solidification  with 
superimposed  predicted  solidification  conditions  for  a  range  of 
withdrawal  rates.  The  simulated  solidification  condition  for  the 
predicted  preferential  withdrawal  rate  is  identified  as  a  grey  circle, 
while  withdrawal  rates  slower  and  faster  than  the  preferred  are 
colored  white  and  black,  respectively. 

A  similar  analysis  was  performed  for  bar  thicknesses  varying 
from  1.3  to  12.7  cm  in  diameter  and  at  furnace  temperatures  from 
1475  to  1600  °C.  The  same  trend  of  maximum  axial  thermal 
gradient  at  an  intermediate  withdrawal  rate  was  observed  (Figure 
9a).  More  specifically,  the  maximum  axial  thermal  gradient  was 
achieved  when  the  position  of  the  solid-liquid  interface  at  the 
exterior  of  the  casting  was  approximately  0.5  to  1.5  cm  and  0.2  to 
0.6  cm  above  the  transition  region  for  the  Bridgman  and  LMC 
processes,  respectively,  independent  of  furnace  temperature  or 
casting  geometry  (Figure  9b). 
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(a)  Withdrawal  Rate  (mm-min ')  Relative  Position  (cm) 

Figure  9:  Predicted  axial  thermal  gradient  for  a  range  of  bar 
thicknesses  cast  at  varying  withdrawal  rates  plotted  as  a  function 
of  (a)  withdrawal  rate  and  (b)  solid-liquid  interface  position 
relative  to  the  furnace  transition  region. 


The  thermal  conditions  associated  with  the  withdrawal  rate  that 
provided  the  maximum  axial  thermal  gradient  for  each  bar 
thickness  or  mold-heater  temperature  mold  configuration  were 
compared  and  superimposed  on  a  solidification  defect  map 
(Figure  10).  Increasing  the  bar  thickness  for  either  the  Bridgman 
or  LMC  process  refined  the  dendritic  structure,  increasing  both 
the  axial  thermal  gradient  and  axial  solidification  rate  without 
proceeding  to  the  nucleated  grain  regime.  For  very  large  bar 
thicknesses  (greater  than  7.5  cm  in  diameter)  produced  via  the 
Bridgman  process,  the  thermal  conditions  are  consistent  with  the 
formation  of  freckles  according  to  the  defect  map.  Increasing  the 
mold-heater  temperature  for  either  the  Bridgman  or  LMC  process 
increased  the  axial  thermal  gradient  and  decreased  the  axial 
solidification  rate.  Based  on  the  predictions,  by  increasing  the 
mold-heater  temperature,  structure  refinement  did  not  occur,  but 
the  propensity  to  form  defects  decreased. 
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Figure  10:  Defect  map  for  directional  solidification  with 
superimposed  predicted  preferred  solidification  conditions  for  a 
range  of  bar  thicknesses  and  mold-heater  temperatures. 

Overall,  maximization  of  the  axial  thermal  gradient  required  a 
balance  of  efficient  heat  extraction  with  minimal  interface 
curvature.  In  addition,  a  coincidence  of  relative  solid-liquid 
interface  position  and  maximum  axial  thermal  gradient  occurred 
for  a  broad  range  of  bar  thicknesses  and  mold-heater 


temperatures.  The  insensitivity  of  the  solid-liquid  interface 
position  to  bar  thickness  and  mold-heater  temperature  suggested 
that  there  is  a  desirable  solidification  position  between  the  hot  and 
cold  zones  of  the  furnace.  Thus,  the  axial  thermal  gradient  was 
utilized  as  the  critical  solidification  metric  and  will  be  the  focus  of 
the  remainder  of  this  paper.  Typically,  the  withdrawal  rate  is  the 
process  parameter  that  is  varied  to  achieve  single-crystal 
solidification.  The  approach  for  each  casting  configuration  herein 
will  be  to  (1)  vary  the  withdrawal  rate,  (2)  maximize  the  thermal 
gradient,  (3)  monitor  the  relative  position  of  the  solid-liquid 
interface  and  (4)  compare  findings  to  classical  defect  criteria. 

Initial  demonstration  of  the  optimization  approach  has  been 
developed  by  comparison  to  experimental  efforts  previously 
published  [17].  Metallographic  analysis  of  castings  from  a  range 
of  withdrawal  rates  demonstrated  that  a  minimum  primary 
dendrite  arm  spacing  (PDAS)  was  achieved  at  an  intermediate 
withdrawal  rate  (Figure  11)  [17].  Superposition  of  the  PDAS 
measurements  and  predicted  axial  thermal  gradients  as  a  function 
of  withdrawal  rate  confirmed  that  the  minimum  PDAS  was 
observed  for  the  experimental  condition  that  provided  the 
predicted  maximum  axial  thermal  gradient,  relative  to  the  other 
experimental  conditions.  Therefore,  not  only  does  the 
optimization  protocol  provide  conditions  that  balance  enhanced 
heat  extraction  and  solidification-front  curvature,  but  also  provide 
optimal  refinement  in  dendrite  structure. 
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Figure  11:  Comparison  of  predicted  axial  thermal  gradient  and 
measured  primary  dendrite  arm  spacing  for  bars  composed  of 
Rene  N5  cast  at  varying  withdrawal  rates  (measurements  adopted 
from  [17]). 


Case  B:  lateral- growth  casting 

Further  analysis  was  conducted  on  the  morphology  of  the 
solidification  structure  and  its  sensitivity  to  the  solid-liquid 
interface  curvature.  This  analysis  (1)  complements  the  previous 
analysis  by  decoupling  the  section  thickness  effect  and  thermal 
conditions  with  respect  to  dendrite-growth  morphology  and  (2) 
highlights  the  importance  of  considering  the  axial  thermal 
gradient  and  axial  solidification  rate  in  the  presence  of  non-axial 
heat  extraction. 


A  previous  analysis  [16]  developed  a  relationship  between 
inclination  angle  and  dendrite-growth  morphology  (Figure  12). 
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More  specifically,  the  onset  of  lateral  growth  occurred  at  an 
inclination  angle  of  25  degrees,  while  the  presence  of  nucleated 
grains  was  observed  at  inclination  angles  greater  than  45  degrees. 
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Figure  12:  Predicted  inclination  angle  of  the  solidification  front 
for  a  range  of  wrap  thicknesses,  withdrawal  rates  and  alloys  with 
superimposed  experimentally  observed  dendrite  morphology  [16]. 


The  thermal  conditions  from  this  analysis  were  compared  to  a 
defect  map  to  assess  whether  there  is  an  underlying  relationship 
with  thermal  gradient  and  solidification  rate,  rather  than 
inclination  angle,  for  the  observed  dendrite-growth  morphologies 
(Figure  13).  The  curvature  of  the  solid-liquid  interface  has  not 
been  previously  considered  with  respect  to  the  solidification 
defect  maps.  The  onset  of  lateral  growth  occurred  at  an 
approximately  constant  G/V  value  of  approximately  3000  °C-s- 
cm'1,  while  grain  nucleation  was  only  observed  at  G/V  values  less 
than  1000  °C-s-cm'1.  These  data  also  indicated  that  lateral  growth 
is  a  precursor  to  the  presence  of  nucleated  grains  as  identified  on 
the  defect  map  (Figure  13).  No  trend  between  G/V  and  dendrite 
morphology  was  observed  when  the  total  thermal  gradient  and 
total  solidification  rate  were  considered. 
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Figure  13:  Defect  map  for  directional  solidification  with 
superimposed  predicted  solidification  conditions  above  the  mold 
wrap  for  a  range  of  wrap  thicknesses  and  withdrawal  rates  to 
achieve  a  broad  range  of  solidification-front  curvature.  The  onset 
of  lateral  growth  and  nucleated  grains  occur  at  a  constant  G/V  in 
the  axial  direction. 


Case  C:  tri-crystal  ring  segment 

In  the  final  analysis,  a  more  complex  geometry  was  simulated  to 
assess  the  validity  of  the  optimization  protocol  to  define  process 
conditions  without  any  experimental  trials.  This  mold 
configuration  was  significantly  different  than  the  rotationally- 
symmetric,  thin- section  castings  considered  thus  far.  A  simulated 
parametric  analysis  determined  the  withdrawal  rate  that  provided 
the  maximum  axial  thermal  gradient,  thus  defining  the  preferred 
process  conditions  for  two  thicknesses  and  two  directional- 
solidification  processes.  Following  the  definition  of  processing 
conditions,  casting  trials  were  conducted  at  the  preferred  process 
conditions  to  evaluate  dendrite  morphology  and  micro  structure 
scale.  Thermocouple  measurements  from  the  casting  trials  were 
also  utilized  to  validate  model  predictions. 

The  predicted  axial  thermal  gradient  at  the  surface  of  the  tri¬ 
crystal  resulted  in  the  same  trend  as  observed  in  the  previous  two 
cases:  a  maximum  axial  thermal  gradient  was  attained  at  an 
intermediate  withdrawal  rate  for  both  thicknesses  and  both 
directional- solidification  processes  (Figure  14).  As  observed 
previously,  the  axial  thermal  gradient  was  less  sensitive  to 
changes  in  withdrawal  rate  for  the  Bridgman  process.  The 
predicted  axial  thermal  gradient  for  the  Bridgman-processed  thick 
tri-crystal  was  near  the  20  °C-cm'1  threshold  for  all  investigated 
withdrawal  rates,  suggesting  that  the  thermal  conditions  may  be 
insufficient  for  single-crystal  growth. 
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Figure  14:  Predicted  axial  thermal  gradient  for  the  thick  and  thin 
variants  of  the  tri-crystal  ring  segment  processed  via  the 
Bridgman  and  LMC  processes,  (a)  as  a  function  of  withdrawal 
rate  and  (b)  as  a  function  of  solid-liquid  interface  position  relative 
to  the  furnace  transition  region. 
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In  order  to  predict  the  dendrite  morphology,  the  predicted 
solidification  conditions  were  superimposed  on  the  classical 
defect  map  (Figure  15).  As  identified  from  the  axial  thermal 
gradient  predictions,  the  Bridgman-processed  thick  tri-crystal  was 
in  the  regime  for  freckle  formation.  The  optimal  process 
condition  for  the  remaining  pairs  of  tri-crystal  configurations  and 
processes  were  in  the  dendrite-growth  regime.  For  each 
configuration  and  directional-solidification  process,  the  path  of 
the  predictions  with  increased  withdrawal  rate  was  consistent  with 
the  5 -bar  cluster  predictions  previously  plotted.  Increased 
refinement  was  predicted  for  the  thin  tri-crystal  casting  as 
compared  to  the  thick,  which  was  consistent  with  the  observed 
trend  of  increased  refinement  with  reduced  cross-sectional 
thickness  due  to  reduced  thermal  mass  and  enhanced  heat 
extraction.  Moreover,  considerable  refinement  was  predicted  for 
LMC-processed  tri-crystals  as  compared  to  the  same  casting 
configurations  processed  via  Bridgman. 
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Figure  15:  Defect  map  for  directional  solidification  with 
superimposed  predicted  solidification  conditions  for  thick  (5.1 
cm)  and  thin  (1.9  cm)  tri-crystal  ring  segments  solidified  via  the 
Bridgman  and  LMC  processes.  The  optimal  solidification 
condition  for  each  pair  of  casting  configuration  and  process  are 
shaded  grey. 


The  modeling-defined  process  conditions  were  utilized  for  a 
casting  trial  for  each  tri-crystal  configuration  for  both  directional- 
solidification  processes.  Three  thermocouples  were  positioned  at 
the  mold-metal  interface  at  the  surface  of  the  casting  at  three 
different  heights  within  the  tri-crystal  segment.  The  predicted 
cooling  rate  was  within  0.05  °C/sec  of  the  cooling  rate  from  the 
thermocouple  measurements  for  each  of  the  four  different  casting 
trials.  PDAS  and  SDAS  measurements  from  similar  locations 
within  each  of  the  tri-crystal  castings  were  compared  to  the 
predictions  of  PDAS  and  SDAS  at  the  same  locations  using  well- 
accepted  dendrite-growth  models  (Table  5)  [20,21].  Predictions 
of  PDAS  and  SDAS  were  within  15  pet  of  the  measurements, 
which  is  considered  to  be  within  the  local  variability  due  to  the 
measurement  techniques  [17]. 


Table  5:  Comparison  of  measured  and  predicted  PDAS  and  SDAS 
for  thick  and  thin  tri-crystal  configurations  processed  via  the 
Bridgman  and  LMC  techniques. 


Casting  Mode 

Section  Thickness  (cm) 

Bridgman 

5.1  1.9 

LMC 

5.1 

1.9 

_ _  _  .  .  Measured 

406 

334 

248 

187 

PDAS  urn 

Predicted 

420 

334 

210 

158 

_ _  _  .  Measured 

86 

72 

43 

35 

SDAS  ( urn) 

Predicted 

89 

72 

49 

41 

Extension  to  Airfoil  Geometries 


With  the  relatively  successful  application  of  the  protocol  to  the  tri¬ 
crystal  ring  segment,  a  hypothetical  airfoil  geometry  was 
considered.  Due  to  the  significant  changes  in  cross-section  along 
the  axis  of  the  blade,  the  optimization  protocol  was  extended  to 
define  a  time-dependent  withdrawal  rate  by  considering  all  casting 
nodes  along  the  blades  axis.  The  time -dependent  optimization 
protocol  has  been  defined  as  follows:  (1)  Conduct  simulations  at  a 
broad  range  of  constant  withdrawal  rates,  (2)  Analyze  the 
simulations  individually  at  all  casting- surface  nodes  for  the 
position  of  the  solid-liquid  interface  relative  to  the  transition 
region  in  the  furnace  as  a  function  of  height  along  the  blade,  (3) 
generate  a  contour  plot  of  relative  solid-liquid  interface  position 
as  a  function  of  axial  position  along  the  blade  and  withdrawal  rate, 
(4)  determine  withdrawal  rate  as  a  function  of  time  by  maintaining 
the  solid-liquid  interface  at  a  specified  relative  position  along  the 
length  of  the  casting. 

This  protocol  was  applied  to  a  CFM56  blade  with  the  following 
results.  The  predicted  relative  position  of  the  solid-liquid 
interface  for  each  incremental  distance  along  the  withdrawal 
direction  was  determined  for  a  constant  withdrawal  rate 
simulation,  Figure  16.  Every  node  at  the  casting  surface  was 
included  in  the  statistical  analysis  (step  2),  but  only  the  average 
value  for  each  0.1  cm  “bin”  was  used  for  the  contour  plot  (step  3). 
An  increase  in  the  relative  position  is  observed  in  transitioning 
from  the  relatively  thick  cross-section  of  the  ramp  to  the  relatively 
thin  airfoil  walls.  A  decrease  in  the  relative  position  is  observed 
as  the  solidification  front  approaches  the  platform  in  the  airfoil 
and  continues  to  decrease  through  the  platform.  The  previous 
analyses  indicated  that  a  maximum  axial  thermal  gradient 
occurred  at  a  relative  position  of  0.5  cm  above  the  floating  baffle 
for  Figure  16.  Thus,  the  normalized  withdrawal  rate  of  1.0 
resulted  in  a  relative  position  that  was  too  high  within  the  ramp 
and  airfoil  regions  of  the  casting,  suggesting  that  an  increase  in 
the  withdrawal  rate  is  required  for  optimal  processing  conditions 
in  these  regions. 
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Figure  16:  Solid-liquid  interface  position  at  the  casting  surface 
relative  to  the  free  surface  of  the  coolant  as  a  function  of  distance 
along  the  casting  for  a  constant  normalized  withdrawal  rate  of  1.0. 


After  the  analysis  of  each  constant-withdrawal-rate  simulation 
was  completed,  the  data  was  combined  onto  a  contour  plot,  Figure 
17.  The  average  value  of  the  relative  position  was  utilized  to 
generate  the  contours.  In  order  to  maintain  a  relative  position  of 
0.5  cm  throughout  the  entire  casting,  a  normalized  withdrawal  rate 
increasing  from  1.2  to  1.5  in  the  ramp  region  and  decreasing  from 
1.5  to  1.1  in  the  airfoil  region  is  required,  which  is  consistent  with 
the  analysis  of  Figure  16,  that  the  normalized  withdrawal  rate  of 
1.0  is  too  slow  to  provide  the  optimal  position. 


Figure  17:  Constant  solid-liquid  interface  position  contours  as  a 
function  of  position  along  the  casting  and  normalized  withdrawal 
rate.  The  solid  black  curve  is  the  optimal  relative  position. 


This  analysis  has  demonstrated  that  it  is  feasible  to  define  a  time- 
dependent  withdrawal  rate  using  only  solidification-modeling 
predictions,  and  that  the  modeling  approach  is  sufficiently  robust 
to  compensate  for  local  variability  based  on  complex  geometry. 
However,  validation  that  the  model-defined  processing  conditions 
achieve  better  solidification  conditions  than  experience-based 
processing  conditions  from  industry  is  ongoing. 


Discussion 

Key  findings  from  this  research  include  (1)  the  development  and 
application  of  the  maximum  axial  thermal  gradient  as  a  thermal 
metric  for  optimization  of  processing  conditions  for  single-crystal 
dendrite  growth,  (2)  an  improved  understanding  of  the  effect  of 
solid-liquid  interface  curvature  on  dendrite-growth  morphology 
and  (3)  the  extension  of  defect-formation  criteria  to  account  for 
curvature  of  the  solid-liquid  interface.  The  implications  of  these 
findings  on  the  processing  of  single-crystal  nickel-base 
superalloys  are  discussed  in  this  section 

The  utilization  of  the  axial  thermal  gradient  for  optimization  of 
the  directional  solidification  conditions  provided  a  balance  of  high 
thermal  gradient  and  reduced  solid-liquid  interface  curvature.  The 
predicted  optimal  solidification  conditions  resided  within  the 
dendrite  growth  regime  of  the  classical  dendrite  morphology  maps 
(Figures  8,  10  and  15).  Furthermore,  a  correspondence  between 
maximum  axial  thermal  gradient  and  solid-liquid  interface 
position  relative  to  the  transition  region  of  the  furnace  indicated 
that  the  appropriate  thermal  metric  was  determined,  since  it  is 
generally  accepted  that  optimal  directional  solidification  occurs 
when  the  solid-liquid  interface  is  in  the  transition  region  of  the 
furnace  (Figures  9b  and  14b).  This  correspondence  reduced  the 
number  of  simulation  required  to  determine  the  optimal 
withdrawal  rate,  since  there  is  a  well-established  relationship 
between  solid-liquid  interface  position  and  withdrawal  rate: 
increasing  withdrawal  rate  decreases  relative  position. 

The  application  of  the  optimization  technique  to  a  complex 
turbine-blade  geometry  required  modification  to  the  optimization 
protocol  in  order  to  allow  for  time-dependent  processing 
conditions.  The  withdrawal  rate  at  each  incremental  distance 
along  the  length  of  the  casting  was  optimized  when  the  average 
relative  position  of  all  of  the  nodes  on  the  casting  surface  is  equal 
to  the  optimal  relative  position.  In  addition,  the  minimum  relative 
position  for  each  incremental  distance  along  the  length  of  the 
blade  was  compared  to  the  metal-coolant  level  to  ensure  that  the 
relative  position  across  the  entire  surface  was  never  below  the  tin 
level  due  to  the  sharp  decay  in  thermal  conditions  below  the 
surface  of  the  metal  coolant.  Micro  structure  scale  and  variability 
measurements  from  casting  trials  using  this  time-dependent 
withdrawal  rate  are  needed  to  confirm  the  validity  of  the 
optimization  approach.  In  addition,  comparison  of  thorough 
experimental  investigations,  such  as  in  Figure  11  and  processing 
experience  from  industry,  such  as  the  dendrite-morphology  maps 
or  casting  geometries  with  optimized  process  conditions,  is  the 
most  effective  means  to  demonstrate  the  value  of  the  optimization 
protocol  on  new  casting  geometries. 

The  experiments  and  simulations  of  the  lateral  growth  mold 
demonstrated  the  need  to  understand  the  role  of  solid-liquid 
interface  curvature  on  dendrite  morphology.  Through  comparison 
of  experimentally  observed  dendrite  morphology  and  predicted 
solidification  conditions,  the  dendrite-morphology  maps  were 
extended  to  account  for  solid-liquid  interface  curvature  through 
utilization  of  the  axial  thermal  gradient  and  axial  solidification 
rate.  The  approximation  that  the  solid-liquid  interface  is  flat  and 
perpendicular  to  the  withdrawal  direction  may  not  be  valid,  even 
for  simple  bar-mold  geometries  as  indicated  by  Figure  12  in 
which  the  solid-liquid  interface  exceeded  20  degrees  of  inclination 
at  the  preferred  withdrawal  rate.  Thus,  future  analyses  should 
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consider  solid-liquid  interface  curvature,  especially  when  utilizing 
FE  solidification  models  capable  of  extracting  the  inclination  of 
the  solidification  front.  An  additional  dendrite-morphology 
regime  was  also  identified  between  misoriented  grain  formation 
and  axial  dendrite  growth. 

Conclusions 

Through  the  iteration  of  solidification  modeling,  the  process 
conditions  can  be  systematically  altered  in  order  to  determine  the 
preferred  process  condition.  The  position  of  the  solid-liquid 
interface  relative  to  the  transition  zone  in  the  furnace  can  be  used 
as  a  metric  to  identify  a  preferred  process  condition.  The  metric  is 
fundamentally  based  on  local  thermal  conditions  at  the  surface  of 
the  casting,  specifically  the  axial  thermal  gradient.  This  metric  (1) 
indirectly  accounts  for  a  balance  between  thermal  gradient  and 
solidification-front  curvature,  (2)  is  independent  of  casting 
geometry  and  furnace  temperature  and  generally  insensitive  to 
alloy,  and  (3)  is  directly  quantifiable  within  solidification - 
modeling  software  for  practical  application. 

Solid-liquid-interface  curvature  affects  the  dendrite  morphology 
and  can  increase  the  propensity  for  misoriented-grain  formation. 
Lateral  growth  has  been  associated  with  increased  inclination  and 
is  a  precursor  to  misoriented-grain  formation. 

Defect-formation  criteria  historically  used  to  define  a  process 
window  were  modified  to  account  for  non-axial  heat  extraction 
and  compared  to  the  solidification  metric  developed  in  this  work. 
The  preferred  process  conditions  defined  by  the  proposed 
solidification  metric  provided  thermal  conditions  that  were  within 
the  recommended  bounds  of  the  classical  solidification  maps. 

A  model  has  been  validated  for  the  directional -solidification 
process  over  a  broad  range  of  thermal  conditions  and  casting 
geometries.  The  optimization  strategy  has  been  verified  on  two 
different  case  histories,  encompassing  a  range  of  casting 
geometries  and  withdrawal  rates  for  two  different  alloys.  Due  to 
its  success,  the  optimization  strategy  was  also  applied  to  a  turbine- 
blade  geometry,  enabling  a  time-dependent  prediction  of  the 
withdrawal  rate  that  provides  for  maximum  axial  thermal  gradient 
by  maintaining  constant  solid-liquid  interface  position  in  the 
transition  region  of  the  furnace. 
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